Quantum interference and supercurrent in multiple-barrier proximity structures 
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We analyze an interplay between the proximity effect and quantum interference of electrons 
in hybrid structures superconductor-normal metal-superconductor which contain several insulating 
barriers. We demonstrate that the dc Josephson current in these structures may change qualitatively 
due to quantum interference of electrons scattered at different interfaces. In junctions with few 
conducting channels mesoscopic fluctuations of the supercurrent are significant and its amplitude 
can be strongly enhanced due to resonant effects. In the many channel limit averaging over the 
scattering phase effectively suppresses interference effects for systems with two insulating barriers. 
In that case a standard quasiclassical approach describing scattering at interfaces by means of 
Zaitsev boundary conditions allows to reproduce the correct results. However, in systems with 
three or more barriers the latter approach fails even in the many channel limit. In such systems 
interference effects remain important in this limit as well. For short junctions these effects result 
in additional suppression of the Josephson critical current indicating the tendency of the system 
towards localization. For relatively long junctions interference effects may - on the contrary - 
enhance the supercurrent with respect to the case of independent barriers. 



I. INTRODUCTION 



In recent years there has been a great deal of activ- 
ity devoted to both experimental and theoretical stud- 
ies of mesoscopic superconducting-normal {SN) hybrid 
structures New important phenomena such as 

anomalous Meissner screening, re-entrant behavior of the 
conductance, nonequilibrium-driven 7r-junction state and 
many others have been discovered and thoroughly inves- 
tigated. In some cases it was found that an interplay 
between the proximity effect and quantum interference 
of electrons in a normal metal play a significant role at 
sufficiently low temperatures. For instance, interference 
of electrons scattered at impurities in the normal layer 
may strongly enhance the Andreev conductance Ga of 
SN systems leading to the so-called zero-bias anomaly 
|^-^| : At low voltages and temperatures Ga turns out to 
depend linearly on the SN interface transmission D in 
contrast to the standard result Ga oc D 2 obtained in the 
absence of interference effects in the iV-layer. 

In this paper we will address a different - although 
somewhat related - problem. We will analyze the dc 
Josephson effect in SNS systems which contain several 
insulating barriers. In this case electrons scattered at dif- 
ferent barriers can interfere inside the junction. We will 
demonstrate that this effect may lead to qualitative mod- 
ifications of the supercurrent across the junction. The 
most pronounced effect of quantum interference is ex- 
pected in SNS systems with few conducting channels. 
This situation can be realized, for instance, if two super- 
conductors are connected via a carbon nanotube [puTl. 
More conventional SNS structures with many conduct- 
ing channels and several insulating barriers are also of 
considerable interest, for instance in relation to possi- 
ble applications, see e.g. Ref. || and further references 
therein. We will demonstrate that for such systems quan- 



tum interference effects are also important provided there 
exist more than two scatterers inside the junction. 

A powerful tool for theoretical studies of mesoscopic 
superconductivity is provided by the quasiclassical for- 
malism of the energy-integrated Eilenberger Green func- 
tions @ (see also Refs. |[||l^Q for a review). The 
Eilenberger equations are, on one hand, much simpler 
than the fully microscopic Gorkov equations and, on the 
other hand, allow to correctly describe the system behav- 
ior at distances much longer as compared to the Fermi 
wavelength 1/fci?- Since typical length scales in super- 
conductors (e.g. the coherence length £o or the London 
penetration depth) are all several orders of magnitude 
greater than I /hp, the quasiclassical approach is usually 
an excellent approximation. 

The quasiclassical equations cannot be applied only 
in the vicinity of inter-metallic interfaces and barriers 
where rapid changes of the system properties (at scales 
comparable to I/Hf) occur. Fortunately, in many cases 
this problem can be circumvented by matching the Eilen- 
berger Green functions on both sides of the interface with 
the aid of the proper boundary conditions. In order to 
derive such conditions it is necessary to go beyond qua- 
siclassics, however under certain assumptions the final 
result can be formulated only in terms of the quasiclas- 
sical Eilenberger propagators. The derivation of these 



boundary conditions was performed by Zaitsev 12 1. Sup- 
plemented by these boundary conditions, the Eilenberger 
quasiclassical formalism was proven to be an extremely 
efficient tool for a quantitative description of numerous 
inhomogeneous and hybrid superconducting structures. 

An important ingredient of the derivation |l2| is the 
assumption that the boundaries are situated sufficiently 
far from each other, so that interference effects emerging 
from scattering at different interfaces can be totally ne- 
glected. Under this assumption one arrives at the nonlin- 
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ear matching conditions involving the third power of the 
quasiclassical propagators. These matching conditions 
are expressed in terms of the interface transparency co- 
efficients for the electrons with different directions of the 
Fermi momenta. It is also essential that Zaitsev bound- 
ary conditions do not depend on scattering phases at the 
interface potentials. 

Although for metallic structures containing one inter- 
face one can indeed disregard interference effects, in sys- 
tems with several boundaries this is in general not any- 
more possible. Hence, the applicability of the nonlinear 
matching conditions [ fl2j to multiple-barrier systems re- 
quires additional analysis. Some authors [fl3| argued that 
the standard quasiclassical approach can break down in 
a multiple-interface geometry due to the problems with 
the normalization of the Eilenberger functions. 

In principle the above problem with boundary condi- 
tions can be avoided within the approach based on the 
Bogolyubov-de Gennes equations [[u]]. However, this ap- 
proach, though frequently successful, may also be techni- 
cally inconvenient in complicated situations, for instance 
because of a necessity to evaluate the energy eigenvalues 
of the system and to perform summation over the energy 
spectrum in the final results. 

It is also possible to formulate an alternative quasiclas- 
sical approach |l5| which allows to avoid the abovemen- 
tioned problems. Without going into details here let us 



just mention that within the technique 1 15 one deals with 
the quasiclassical spinor amplitude u, w-functions which 
depend on one coordinate and one time only and obey 
linear first order equations. The Eilenberger Green func- 
tions are expressed via two linearly-independent solutions 
of these equations in a way that both the Eilenberger 
equations and the normalization conditions are automat- 
ically satisfied. The formalism (TH) - just as the Eilen- 
berger one - can be formulated both within the Mat- 
subara and the Keldysh techniques and thus is suitable 
both in equilibrium and non-equilibrium situations (see, 
e.g. Ref. r|) in various superconducting structures. It 
is also important that very simple linear boundary condi- 
tions for the quasiclassical amplitudes can be formulated 
at each of the interfaces where electron scattering takes 
place. The number of interfaces in the system is not re- 
stricted and the interference effects are properly taken 
care of. Thus it is possible to take advantage of the qua- 
siclassical approximation and at the same time to for- 
mulate general and simple boundary conditions without 
making additional assumptions employed in Ref. f^ j. 

Similar ideas have recently been put forward by She- 
lankov and Ozana |l7j]. These authors also used linear 
matching conditions (obtained by means of the scattering 
matrix approach) for the "wave functions" which factor- 
ize the two-point Green functions. The next step JlT[ ] was 
to construct quasiclassical one-point Green functions and 
formulate nonlinear boundary conditions for such func- 
tions which would now adequately include information 



about scattering on arbitrary number of "knots" . Linear 
boundary conditions were also used by Brinkman and 
Golubov [jl8| in a calculation related to ours (see below). 

In this paper, following Refs. |lj,[l7],[l8| , we will use 
simple linear boundary conditions in order to match the 
quasiclassical amplitude functions at interfaces. How- 
ever, unlike in Ref. Jl7t , we will avoid reformulating these 
boundary conditions as nonlinear ones for the Eilenberger 
Green functions. Rather we will directly express the 
two-point Green functions and the expectation value of 
the current operator in terms of the quasiclassical am- 
plitudes. We will then apply our method to the calcu- 
lation of dc Josephson currents in hybrid SINI'S and 
SINI' NI" S structures in the clean limit and for arbi- 
trary interface transmission coefficients 

The interference of the scattering events at different 
interfaces manifests itself in the expressions containing 
scattering phases <f> a t the interface potentials. For the 
systems with two barriers (in our case S/iV/'S'-systems) 
with many transmission channels the summation over 
their contributions is equivalent to effective averaging 
over <f>. In this limit one can demonstrate that after 
such averaging our result is equivalent one obtained from 
the Eilenberger equations supplemented by the Zaitsev 
boundary conditions. However, in the case of more than 
two barriers (i.e. for SINI' NI" S junctions) the depen- 
dence on the scattering phases turns out to be much more 
essential. In this situation the approach employing Za- 
itsev boundary conditions turns out to fail also in the 
many channel limit where quantum interference effects 
survive even after averaging over the scattering phases. 

The paper is organized as follows. Our quasiclassical 
approach is outlined in Sec. II. In Sec. Ill we apply 
this approach for the analysis of the dc Josephson effect 
in SINI' S structures with arbitrary interface transmis- 
sions. The Josephson current across SINI' NI" S struc- 
tures is evaluated in Sec. IV. In Sec. V we present a brief 
discussion and summary of our results. Some technical 
details of our calculation are relegated to Appendices. 



II. GENERAL METHOD 

A. Quasiclassical approximation 

The starting point of our analysis are the microscopic 
Gor'kov equations |19|| . In what follows we will assume 
that our system is uniform along the directions parallel 
to the interfaces (coordinates y and z). Performing the 
Fourier transformation of the normal G and anomalous 
F + Green function with respect to these coordinates 



G Un (r,r') 



d 2 k 



(2tt 



\G Un {x,x' ^e^^ 



we express the Gor'kov equations in the following stan- 
dard form 
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iui n — H A(x) 
A* (a;) iuj„ + H c 



G Un (x,x',k\\) \ _ 
F£(x,x',k\\ 



5(x-x') 




(1) 



Here uj n = (2n + 1)tvT is the Matsubara frequency, 
and A(x) is the superconducting order parameter. The 
Hamiltonian H in Eq.fffl) reads 



H 



2m dx 



k 2 

2m 



e F + V(x). 



(2) 



Here fey = fcy — -A||(x), cf is Fermi energy, the term 
V(x) accounts for the external potentials (including the 
boundary potential). The Hamiltonian H c is obtained 
from H (Q) by inverting the sign of the electron charge e. 
The above Hamiltonians can also include the self-energy 
terms which, however, will not be considered below. 

The quasiclassical approximation makes it possible to 
conveniently separate fast oscillations of the Green func- 
tions due to the factor exp(±ik x x) from the envelope of 
these functions changing at much longer scales as com- 
pared to the atomic ones. Making use of this approxi- 
mation for two-component vector Tp ± (x) exp(±ik x x) we 
obtain 



iuj n — H 
A*(x) 



A(x)^ 

tCJn + H r 



ip±(x)e 



iuj n - HI 
A*(x) 



A(x) 



H a ±c 



<P±(x), ( 3 ) 



where we defined k x 



fen and 



fl± = Tiv x d x - -A\\(x)v\\ 



2mc 2 



A l[ 2 (x) + V(x). (4) 



Here v x — k x /m, V(x) represents a slowly varying part 
of the potential which does not include fast variations 
which may occur at metallic interfaces. The latter will 
be accounted for by the boundary conditions to be formu- 
lated below. But first let us briefly describe the general 
structure of the Green functions obeying eq. (Q) . 

B. Construction of the Green functions 



Consider the equation 



lUJ n 



HI 



A*(x) 



A(x) 
iuj n + H^ c 



= 0. 



(5) 



There exist two linearly independent solutions Tp + of eq. 
(H). One such solution (denoted below by Tp +1 ) docs 
not diverge at x — > +oo, the other solution ~ip +2 is well- 
behaved at x — > — oo. Similarly, two linearly independent 
solutions Tp_^ 2 do not diverge respectively at x — > — oo 
and x — * +oo. 



A particular solution of the Gor'kov equations (|]) can 
now be sought in the following form 



Tp_ 2 {x)g 2 {x')e^ kAx - x ' ) ifx>x' 



(6) 



and 



(?££:^!))=^- i( ^^ )e "^ (a: " x ' )+ 

Tp +2 {x)f 2 {x')e lk ^ x - x '^ ifx<x'. (7) 

These functions satisfy Gor'kov equations at x ^ x'. The 
functions fi^(x) and g\ t2 (x) are determined with the aid 
of the continuity condition for the Green functions at 
x = x' and the condition resulting from the integration 
of 8{x — x 1 ) in eq. ([!]). As a result we arrive at the linear 
equations 



t P+i{x)gi(x) + tp_ 2 (x)g 2 (x) = 
T P+x{x)gx{x)-lp_ 2 {x)g2{x) 



V-x{x)fx{x) - v3 +2 (x)/ 2 (x) 



(8) 



which can be trivially resolved. 

For a homogeneous superconductor in the absence of 
the magnetic field this procedure allows to immediately 
recover the well known result 

G Un (x,xT) = - Ml ' +7 *) (e lks \ x " x '\ - 7 2 e -ifeSlx-x'h ( 

F+Jx,x') = v £* 3) ( e «*l*-*'l + e -<*Sk-*'l) , 

where x 1S the phase of the pairing potential, k$ = 
k x + iVL n /v x , Q„ = v/|A| 2 + and 7 = ■ Here 

for convenience we set ui n > 0. 

In a non-homogeneous situation a general solution of 
the Gor'kov equations takes the form 



G Un (x,x')\ _ (G Un {x,x') 
F+(x,x') ) \ F+ (x,x') 



part 



[h(x')ip +1 {x) + l 2 {x')ip +2 {x)]e 1 ' 
[l 3 (x')if^(x) + h{x')lp_ 2 {x)]e-^ x . 



(9) 



For systems which consist of several metallic layers the 
particular solution is obtained with the aid of the pro- 
cedure outlined above provided both coordinates x and 
x' belong to the same layer. Should x and x' belong to 
different layers, the particular solution is zero because in 
that case the (5-function in eq. (|l]) fails. The functions 
2, 3,4(2/) in each layer should be derived from the proper 
boundary conditions which we will now specify. 



3 



C. Boundary conditions 

In what follows we shall assume interfaces to be non- 
magnetic. In this case matching of the wave functions 
on the left and on the right side of a potential bar- 
rier, respectively A\ exp(ikx x x) + B\ exp(— ik% x x) and 
Ai exp(ik2xx) + B2 exp{—ik,2xX), is performed in a stan- 
dard way (see e.g |20f): 



A 2 = aAx+PBx, B 2 = f3*A 1 +a*B 1 , 

|2 I 01 2 ^lx 



or 



k 2x 



(10) 



The reflection and transmission coefficients are given by 



R 



D = l-R 



k2x\a\ 



(11) 



For the sake of simplicity below we shall set k\ x = k 2x . 
Since typical energies of interest, such as A and typi- 
cal Matsubara frequencies, are all much smaller than the 
magnitude of the interface potentials, the relationships 
( |l0| ) can be directly applied to the two-element columns 
in eq. (|^). In this way we uniquely determine the Green 
functions of our problem. 

For an illustration let us consider a metallic layer 
with the left and right boundaries located respectively 
at x — di and x — d 2 . We will also choose the argument 
x' inside this layer. As it was already explained, the par- 
ticular solution of the Gor'kov equations for x < d\ or 
x > c?2 is G(x, x') = F + (x, x') = 0, while it has the form 
©i (@) if the coordinate x belongs to this layer. Thus at 
the left boundary we get 



lp +2 (d 1 )f2(x')e- ik * x ' +h(x')<? +1 (d 1 ) 

ai [iH^liidi) +ii(x')^ 2 (d 1 )j + 
A [z 3 V)^(di) + z 4 V)^- 2 (di)] 

and 



h(x')ip +2 {di) 



(12) 



<P- 1 (di)f 1 (x')e ik '*' + I 3 (aO?U(di) + h(x'p- 2 (di) = 

PI [Jf(a0^i(di) + #(aO?^ 2 (di)] + 

a\ ^(x')^-i(di) + ZfV)P- 2 (di)] ■ (13) 

The superscript L labels the solutions in the layer lo- 
cated at a; < d\. The above boundary conditions pro- 
vide four linear equations for the functions l(x') with the 
source term f\, 2 (x') exp(±ik x x'). Similarly, with the aid 
of eq. (^) four boundary conditions at the right boundary 
x = d 2 can be established. Analogous procedure should 
be applied to other interfaces. 



III. JOSEPHSON CURRENT IN SINI'S 
JUNCTIONS 

We shall consider SNS junctions composed of clean 
superconducting (S) and normal (N) metals. We will 



assume that a thin insulating layer (/) can be present 
at both SN interfaces which, therefore, will be charac- 
terized by arbitrary transparencies ranging from zero to 
one. Specular reflection at both interfaces will be as- 
sumed. We also assume that between interfaces elec- 
trons propagate ballistically and no electron-electron or 
electron-phonon interactions are present in the normal 
metal. For simplicity we will restrict our attention to the 
case of identical superconducting electrodes with singlet 
isotropic pairing. Furthermore, we shall neglect possi- 
ble suppression of the superconducting order parameter 
A in the electrodes close to the SN interface. This is a 
standard approximation which is well justified in a large 
number of cases. The phase of the order parameter is set 
to be — \/2 in the left electrode and %/2 in the right one. 
The thickness of the normal layer is denoted by d. 

In order to evaluate the dc Josephson current across 
this structure we shall follow the quasiclassical approach 
described in the previous section. Technical details of our 
calculation are presented in Appendix A. As a result we 
arrive at the expression for the two point Green function 
in the normal layer. After that the current density can 
be calculated from the standard formula 



m 



(14) 



On can also rewrite the current in the form J = J+(x) — 
J+(~x) where J+ is defined by eq. (ph with positive 



Matsubara frequencies u> n > 0. Using (A9) and omitting 
terms oscillating at atomic distances we obtain 



J+ = 2ieT ( 

Lo n >a-'\ k \\\< k '- 



J^{V 1 -U2), (15) 



or explicitly 



J = 4eTsinx 



uj n >Q 



k x dk x sinx 
2n cos x + W ' 



(16) 



where we defined 



ma 



w = 



DiD 2 A 2 
R 2 )+ujID 1 D2 



cos(2k x d 



+ 



D!D 2 A 2 
2(1 - RxRa) n n u n 



D 1 D 2 



A 2 



cosh 



sinh 



2u n d 

v x 
2uj n d 

v x 



(17) 



Here 2k x d + (f> is the phase of the product a 2 /3 2 a*/3*. 
Eqs. (|l^), ( |T7| ) provide a general expression for the dc 
Josephson current in SINI'S structures valid for arbi- 
trary transmissions D± and D2 ranging from zero to one. 
This expression is the central result of this section. 
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We also note that the integral over k x in eq. ( |l6| ) can 
be rewritten as a sum over independent conducting chan- 
nels 



A_ 

2- 



k F N 

k x dk x (...) -» 



(18) 



where A is the junction cross section. In this case D\2 
and i?i,2 may also depend on m. This would correspond 
to different transmissions for different conducting chan- 
nels. 

Finally let us point out that in the limit of symmetrical 
low transparent barriers D\ = D2 <C 1 the problem was 
recently studied by Brinkman and Golubov | fL8| . In the 
corresponding limit their result (eq. (8) of Ref. |fl8|) i s 
similar - although not fully equivalent - to our eqs. (16), 

(0)- 



A. Junctions with few conducting channels 

Let us first analyze the above result for the case of 
one conducting channel N = 1. We observe that the first 
term in eq. ( |l7| ) contains cos(2k x d+4>) which oscillates at 
distances of the order of the Fermi wavelength. Provided 
at least one of the barriers is highly transparent and/or 
(for sufficiently long junctions d > £0) the temperature is 
high T vp/d this oscillating term is unimportant and 
can be neglected. However, at lower transmissions of 
both barriers and for relatively short junctions d < vp/T 
this term turns out to be of the same order as the other 
contributions to W (0). In this case the supercurrent 
is sensitive to the exact positions of the discrete energy 
levels inside the junction which can in turn vary consid- 
erably if d changes at the atomic scales ~ \ jkp. Hence, 
one can expect sufficiently strong sample-to-sample fluc- 
tuations of the Josephson current even for junctions with 
nearly identical parameters. 

Let us first consider the limit of relatively short 
SI NTS junctions in which case we obtain 



e A T sin \ 



V 



tanh 



X>A 
~2T 



(19) 



where we defined 



P = ^/l-Tsin%/2) (20) 

and an effective normal transmission of the junction 

D X D 2 



T = 



1 + R1R2 + 2ViOTcos(2fc x <2 - 



(21) 



Eq. ( |l9| ) has exactly the same functional form as the re- 
sult derived by Haberkorn et al. ]2l[ for SIS junctions 
with an arbitrary transmission of the insulating barrier. 



This result is recovered from our eqs. fll9|), ( |2l| ) if we as- 
sume e.g. D\ <C D2 in which case the total transmission 
((2l]) reduces to T ~ D x . 

As we have already discussed the total transmission T 
and, hence, the Josephson current fluctuate depending on 
the exact position of the bound states inside the junction. 
The resonant transmission is achieved for 2k x d+<p = ±7r, 
in which case we get 



T — 

1 ros 



D X D 2 



(22) 



(1 - x/iOT) 2 ' 

This equation demonstrates that for symmetric junctions 
Di = D 2 at resonance the Josephson current does not 
depend on the barrier transmission at all. In this case 
T ICS = 1 and our result (|l9| ) coincides with the for- 
mula derived by Kulik and Omel'yanchuk |22| for bal- 
listic constrictions. In the limit of low transmissions 
Dx,2 < 1 we recover the standard Breit-Wigner formula 
T Ies = AD 1 D 2 /(D 1 + D2) 2 and reproduce the result ob- 
tained by Glazman and Matveev |23) for the problem of 
resonant tunneling through a single Anderson impurity 
between two superconductors. 

Note that our results (|l9|-pl"|) also support the conclu- 
sion reached by Beenakker p4| that the Josephson cur- 
rent across sufficiently short junctions has a universal 
form and depends only on the total scattering matrix of 
the weak link which can be evaluated in the normal state. 
Although this conclusion is certainly correct in the limit 
d — > 0, its applicability range depends significantly on 
the physical nature of the scattering region. From eqs. 
(|l6|), ( |i~7| ) we observe that the result fll9|) , (pQ) applies 
at d <C £0 not very close to the resonance. On the other 
hand, at resonance the above result is valid only under 
a more stringent condition d <C £oZ?max, where we define 
-Dmax =msx(Di,D 2 ). 

Now let us briefly analyze the opposite limit of suf- 
ficiently long junctions d 3> Co- Here we will restrict 
ourselves to the most interesting case T = 0. From eqs. 
@, @ we obtain 



ev x sin x 
irdzi 



arctan yj z 2 /z\ 



\Jz2lz1 



(23) 



Zl =cos 2 ( X /2) + 
z 2 =sin 2 (x/2) + 



1 



DiD 2 
1 



R+ + 2\/R 1 R 2 cos(2k x d + 4>) 



R+ - 2\/R 1 R 2 cos(2k x d + <p) \ , 



D X D 2 

where i?+ = R± + R 2 . For a fully transparent channel 
D% = D 2 = 1 the above expression reduces to the well 
known Ishii-Kulik result |23j2C 



/ = 



ird ' 



-7T < X < A", 



(24) 



whereas if one transmission is small Oi < 1 and D 2 w 1 
we reproduce the result |27| 

ev x Di sinx 



2d 



(25) 
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Provided the transmissions of both ./VS'-interfaces are low 
-Di,2 C 1 we obtain in the off-resonant region 



I=-^D x D 2 smxf[2k x d + <t>], 
47ra 

where T[x] is a 27r-periodic function defined as 

x 



T[x] 



-7T < X < 7T. 



(26) 



(27) 



In the vicinity of the resonance ||2A x d + <p\ — 7r| < -D max 
the above result does not hold anymore. Exactly at res- 
onance 2k x d + (j) = ±7r we get 



from all the channels. Although all these contributions 
have the same form, they are in general not equal be- 
cause the phase factors 2k x d + (j> change randomly for 
different channels. Accordingly, mesoscopic fluctuations 
of the supercurrent should become smaller with increas- 
ing number of channels and eventually disappear in the 
limit of large N. In the latter limit the Josephson cur- 
rent is obtained by averaging over all values of the phase 
2k x d+(j). The corresponding results are presented below. 



ev x y/DiD 2 sinx 




(28) 



For a symmetric junction D± 2 = D this formula yields 



ev x Dsm(x/2) 
2d 



-TT < X < TT, 



(29) 



while in a strongly asymmetric case D\ <C D 2 we again 
arrive at the expression (p5f). This implies that at reso- 
nance the barrier with higher transmission D 2 becomes 
effectively transparent even if D 2 <C 1. We conclude that 
for £>i,2 <C 1 the maximum Josephson current is pro- 
portional to the product of transmissions D\D 2 off reso- 
nance, whereas exactly at resonance it is proportional to 
the lowest of two transmissions D\ or D 2 . 

We observe that both for short and long SINI'S junc- 
tions interference effects may enhance the Josephson ef- 
fect or partially suppress it depending on the exact po- 
sitions of the bound states inside the junction. We also 
note that in order to evaluate the supercurrent across 
SINI'S junctions it is in general not sufficient to de- 
rive the transmission probability for the corresponding 
NINI'N structure. Although the normal transmission 
of the above structure is given by eq. for all values of 
d, the correct expression for the Josephson current can be 
recovered by combining eq. (21) with the results [ pl|j24| 
in the limit of short junctions d <C -D£o only. In this case 
one can neglect suppression of the anomalous Green func- 
tions inside the normal layer and, hence, the information 
about the normal transmission turns out to be sufficient. 
On the contrary, for longer junctions the decay of Cooper 
pair amplitudes inside the A^-layer cannot be anymore 
disregarded. In this case the supercurrent will deviate 
from the form ( J19|) even though the normal transmission 
of the junction (|2 1|) will remain unchanged. This devia- 
tion becomes particularly pronounced for long junctions, 
i.e. for d ^> £ out of resonance and for d 3> D^ at 
resonance. 

The above analysis can trivially be generalized to the 
case of an arbitrary number of independent conducting 
channels inside the junction N > 1. In that case the su- 
percurrent is simply given by the sum of the contributions 



B. Many channel limit 



Averaging over the phase factors 2k x d + <j> is effectively 
performed by integrating over directions of the electron 
momentum in eq. (|lrj). Since the term in the expression 
for W (|l7j ) which contains cos(2k x d + 4>) oscillates very 
rapidly with changing k x , averaging can be performed by 
first integrating the current ( |l6|) over the phase 2k x d + <fi 
and then integrating the result over k x . We obtain 



— ek F T sin x 

w„>0 



1 hj^jfi) 



(30) 



Here and below we define [i = k x /k,F, ti j2 (ii) 
D lt2 (fj,)/(R lt2 (fx) + 1), i± = ti ± t 2 and 



Q 



ht 2 cosx+ (l + (M 2 + l)^f)cosh- 



uj n fl n . 2uj n d 
+1+ — r"o Sinn 



n 2 



A 2 



flV F 



ptv F 

(1 -*?)(i 



(31) 



The above equations fully determine the Josephson cur- 
rent in SINI'S junctions in the many channel limit and 
at arbitrary transmissions of specularly reflecting SN- 
interfaces. 

Let us make use of this result in order to perform a di- 
rect comparison between our analysis and the approach 
based on the Eilcnberger equations supplemented by Za- 
itsev boundary conditions. The corresponding calcula- 
tion within the latter approach is performed in Appendix 
B. It is interesting to observe that for SINI'S junctions 
this calculation yields exactly the same result (|3^), ( |3l| ) 
as obtained within our calculation after averaging over 
the scattering phase 2k x d + 4>. 



G 




Fig. 1. The Josephson current density (j33|) normalized 
by ekpVF /6ir 2 d is plotted as a function of the phase differ- 
ence X- Here we assumed that the boundary is described 
by an effective potential Uo8(x ± d/2) in which case one has 
t = Vx/(2Uq + The dependence J(x) was evaluated for 
2U§/v F = 10" 4 , 0.03, 0.1, 0.3, 0.6 (from top to bottom). 

This observation allows to make an important con- 
clusion concerning the applicability of the quasiclassical 
analysis employing Zaitsev boundary conditions for the 
Eilcnbcrger propagators. The exact result for the Joseph- 
son current in SINI'S systems, eqs. (|l6|), jl7|), cannot 
be recovered within the latter approach because it es- 
sentially ignores interference effects arising from electron 
scattering on two insulating barriers. At the same time, 
in the limit of many conducting channels the scattering 
phase is effectively averaged out. In this limit Zaitsev 
boundary conditions turn out to correctly describe the 
supercurrent. It is also important to emphasize that the 
latter conclusion applies for the systems with not more 
than two barriers. Below we will analyze the supercur- 
rent in SNS structures with three insulating barriers and 
will show that the approach based on Zaitsev boundary 
conditions fails to provide correct results even in the limit 
of many conducting channels. 

But first let us present several limiting expressions for 
the sake of completeness. We start from the limit of a 
sufficiently thick junction rf>(o and consider T = 0. In 
this case we find 



J 



e *^jf * f WWEAFfo h), (32) 



where F(<p, h) = Jq(1 — h sin 2 9)~ 1 / 2 d9 is the incomplete 
elliptic integral, ip = arcsin(l/ v / /7), h = fx + — /1/2 
and 

1 

h = 



1 - ht 2 cos 2 f - i [1 - ht 2 - >/(l-t?)(l-^) 



h = 



1 - ht 2 sin 2 § - i [1 - M2 - V(l-*i)(l-*2)J 



For an SINS junction (t 2 — 1) the above result yields 
(cf. Ref. H) 

ekpVF siny f 1 duu 2 ti ll — 1± cosy 

J = — , zz / — arctan - ' 



dn 2 



yjl - t\ COS 2 % 



1 + 1\ cosx 



The expression (|32|) also simplifies in the case of a sym- 
metric junction t\ = t% 



^ ekjpVF f 1 pjJi 2 djjL 



K 2 d J Q + p- 



■Fly 



1 



1 + ? 



(33) 



where 



y = arccos [t cos(x/2)] , p(fi, x) 



t 2 sin x 



2VT _7 < 2 



The current density J ( p3|) is plotted in Fig. 1 as a func- 
tion of the Josephson phase x f°r several values of the 
barrier transmission. Note that in the case of small inter- 
face transparencies the limit T — > is effectively achieved 
at temperatures much lower than tvF/d. 

Let us now proceed to the case of small transparencies 
of both interfaces t\ t2 <C 1. In this limit the expression 
(|3l|) takes the form 



ft 4 

Q= 
^ A 4 



. 2u n d w„ 
smh h t+ — 



%H^XI (34) 



where 

P(/i,x)=t+(M)cos 2 (x/2)+t 2 _(/i)sin 2 ( X /2). (35) 

As we have already pointed out, the above result is not 
identical to one presented in eq. (13) of Ref. |1| (see 
also B). However, it is easy to see that this difference 
does not affect the final expression for the current in two 
important limits of short (d <C t£o) and long [d 3> i£o) 
junctions. Only in the intermediate case d ~ ££o some 
deviations between our results and those of Ref. [jOI are 
observed. This is demonstrated in Fig. 2. 

The case of short junctions d <C t£o was already studied 
in Ref. [Gl|. Therefore here we only present the asymp- 
totic expression for the current at d 3> ££o 



ekpVF smx 
2^d 



dpp 2 tit 2 ln(ei/e 2 ), 



(36) 



where t\ = minj/iw^/d, A}, e 2 = (ivp/i^dvV) for T <C 
tvp/d and e 2 ~ T for tvp/d <C T <C ei. The accuracy 
of the above formula is in general logarithmic, and it be- 
comes next to logarithmic in the limits d <C £0 or d £o- 
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Fig. 2. The maximum Josephson current J c of a sym- 
metric SIN I S-junction with low transparency. The angular 
dependence of transparency is taken to be t(/i) — tofJ 2 . The 
normal-state resistance of the structure is denoted by Rn- 
The solid line represents eqs. ( jjo| ) , ( ^) of our paper and the 
dashed line represents eq. (13) of Ref. [18]. 



IV. JOSEPHSON CURRENT IN SINI'NI"S 
JUNCTIONS 

Let us now consider SNS structure with three insu- 
lating barriers. As before, two of them are located at 
SN interfaces, and the third barrier is inside the N -layer 
at a distance d\ and d 2 respectively from the left and 
right SN interfaces. The transmission and reflection co- 
efficients of this intermediate barrier are denoted as Dq 
and Rq = 1 — Dq, whereas the left and the right barri- 
ers are characterized respectively by D\ = X — Ri and 
D 2 =l-R 2 - 

The supercurrent is evaluated along the same lines as 
it was done in Section 2 for the case of two barriers. A 
straightforward (although somewhat lengthy and cum- 
bersome) procedure yields the final result for the Joseph- 
son current which is again expressed by eq. (|l^), where 
the function W now takes the form 



W = W+ + W- + W 12 



(37) 



These three contributions to the VF-function depend re- 
spectively on the sum of thicknesses d\ and d 2 , their dif- 
ference and on these two values separately. We find 



W + ^^^^cos[2k x (d 1+ d: 



'i + *) + 01 + hi + ^ff 1 ^ sinh ^ 

n l(l+R 1 )(l+R 2 )+ulD 1 D2 , 2^„(di+d 2 ) . 

DqD\D 2 A 2 L v x ' 

"doDiD^ cos[2fc x (di 



+d 2 ) 



vv - — DoDiD-2 A 



El. 

D 



- 1 , 2R (D 1 -D 2 ) n n uj n • 1 2u„(d 1 -d 2 ) 

<h) + <Pi - <p 2 \ H d„d 1 d 2 " smn 

2fll 



+ ^^cos[2k x d 2 + ^ 2 ] f< 



"l+^^tflal C o S h ^(d!-d a ) 

A 2 D\D 2 v x 

n 



A 2 v x 



D 2 A 2 LUM1 

. sinh 2^k + i±Ri^| C o S h 5i 

v x Di A 2 



(38) 
(39) 

(40) 



Here we introduced two phases 

2k x d 1 +4>i = axgal/3 al/3l, (41) 
2k x d 2 + <t> 2 = argao/35«2/32, 

related to the corresponding elements of the scattering 
matrices for all three barriers. In contrast to the case of 
two barriers these phases cannot be simultaneously re- 
moved by shifting k x . The expression for the Josephson 
current in SIN I' S junctions derived in the previous sec- 
tion can easily be recovered if we set Dq = 1 — Rq = 1. 
By setting D\ t2 = 1 — R\% = in the above equations 
we arrive at the result for the supercurrent in SNINS 
systems derived in Ref. [E^J . 

A. One channel limit 

Let us first analyze the above general result in the limit 
of one conducting channel. In the limit of short junctions 



d <C Co-Dmax we again reproduce the result (|l9| ) where the 
total effective transmission of the normal structure with 
three barriers takes the form 



1 + 41*0*2 +C(V?1,2,*0,1,2)' 

where 

C = COS^^l - tl){\ -t\)+ COS tp 2 y/(l - i§)(l - t\) 



+ (coster cos</?2 — to sin^i sin ^2) y (1 — t\)(l — 

(43) 

Here we define to, 1,2 = -Do, 1,2/(1 + i?o,i,2) and ipi j2 — 
2k x d\^ 2 + 0i,2- For later purposes let us also perform av- 
eraging of this transmission over the phases ipi and ip 2 . 
We obtain 
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(T) = 



2t\tot 2 



(44) 



In particular, in the case of similar barriers with small 
transparencies -Do, 1,2 ~ -D <C 1 the average normal trans- 
mission of our structure is (T) ~ D 3 ! 2 . Suppression of 
the average transmission below the value ~ D is a result 
of destructive interference and indicates the tendency of 
the system towards localization. Eq. ( fl4| ) follows from an 
explicit integration, but it can also be understood in sim- 
ple terms. Consider the square < tp 1 < 2n,0 < <p 2 < 
2t:. The main contribution to the average transmission 
comes from the resonant region T ~ 1. In the symmetric 

case to, 1,2 = t <C 1 this resonance occurs approximately 
r 1 2 

along the lines y(l + cos yi)(l + cos ^2) — i ~ i 3 
in two quadrants ipi,ip% < tt and ipi,if2 > 7r. In 
other words, the resonant region is represented by two 
hyperbola- like curves with characteristic widths ~ D 3 / 2 . 
This dependence of the average transmission is recovered 
from the exact result (B). 

Let us now proceed to the limit of a long junction 
di,2 S> £0 and T = 0. In the off-resonant region we 
find 



J = 3—5 o(Vi,2, oa/ffli)) 

8irdi 



where 

B = r * . 

Jo [coshx + cosy>i] [cosh(d2X / di) + cosy^] 
Evaluating this integral for d\ = d 2 we get 
ev x D 1 D D 2 sinx T[<pi] - T[y> 2 ] 



J = 



8?rdi 



COS </?2 — COS Ifi 



(45) 



(46) 



(47) 



This expression diverges at resonance (i.e. at <pi ~ 7r or 
</?2 — 7r) where it becomes inapplicable. In the resonant 
region ip 2 — tt we obtain 



J = 



ev x y/D 1 D D 2 sinx 



4^2(1 + cos ^ )(T~ 1 - sin 2 ( x /2)) 



(48) 



B. Many channel junctions 

As it was already discussed, in the many channel limit 
it is appropriate to average the current over the scattering 
phases. If the widths d\ and d 2 fluctuate independently 
on the atomic scale, averaging over tpi and tp 2 can also be 
performed independently. If di and d 2 do not change on 
the atomic scale but are incommensurate, independent 
averaging over the two phases can be performed as well. 
The situation is different only for strictly commensurate 
d\ and d 2 in which case no independent averaging can be 
fulfilled. 



Let us first briefly discuss the latter situation of com- 
mensurate A^-layers. For simplicity we assume d\ = d 2 , 
consider a symmetric situation D\ — D 2 =D<1 and 
set the transparency of the intermediate interface to be 
Do S> D 2 . We will only present the result for the case 
of short junctions d <C £o-Dmax- We observe that the de- 
nominator in eq. (|l6|), (37) has a resonant structure as a 
function of ipi + ip 2 . Integrating near the resonances we 
obtain 



1 f 1 

J = — ekpTsinx / dfifiD(fi) Y] 
2lT Jo ,.fV, 



where 



K= (1 + l/JoTsin 2 , 



1 - 



%n 2 ny /n(x^) 



A 2 sin 2 ( X /2) 



^(l + l/3o| 



sin 



(49) 



(50) 



and |/?o| 2 (^) = [1 — Do(fi)]/Do(n). An interesting fea- 
ture of the expression ( p0|) is a dependence of the crit- 
ical Josephson current on the scattering phase <fi — 
(4>i — 4> 2 )/2. For instance, provided \(3q\ is large (the 
transparency of the intermediate layer is small), the crit- 
ical current can vary from ~ ekp\A\D to ~ ekp\A\DDo 
with (f> varying from to ir/2. One should also bear in 
mind that <p may depend on fi. However, since the main 
contribution to the supercurrent comes predominantly 
from the electrons with the momenta perpendicular to 
the interfaces, we can estimate the current with (j) cor- 
responding to the forward direction. If Dq <C D 2 and 
sin 2 4>~^> D the current is given by eq. (|o|) . For = 
(or (j) — 71") a different expression follows 



1 

2~n' 



J = —ekp sin %A tanh 



A 

2T.,n 



(51) 



Now let us consider a more realistic situation of in- 
commensurate d\ and d 2 which allows for independent 
averaging over the scattering phases ipi and ip 2 . Tech- 
nically this procedure amounts to evaluating the inte- 
gral of the expression l/[t + cos x cos (Ax)] from x = 
to some large value x — L. At A = 1 the result of this 
integration is L/y/t(l + 1). However, if A is irrational, 
the integral approaches the value 2LK(l/t 2 )/irt, where 
K(h) = F(ir/2,h) is the complete elliptic integral. This 
simple example illustrates our averaging procedure over 
two independent phases x and Xx. 

Let us assume that the transparencies of all three in- 
terfaces are small as compared to one. After averaging 
over ipi one arrives the expression which has a resonant 
dependence on ip 2 near tp 2 — 7T. Expanding in powers of 
near this resonance with 6<p 2 — tp 2 — tt and keeping the 
terms proportional to Sip 2 and 5ip 2 we find 



cos x+W + +W-+W 12 



A 

21}- 



28tp 2 



D0D1D2 



2A 2 sin 2 (x/2) 

nn 



S V>2 

Dim 
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Then evaluating the integral over i5</?2 we derive the final 
expression for the current 

J = ^fD cS smxT 



A 2 

w„>0 n 



rA 2 sin^( X /2) 



fi2 



where we define the effective transmission 



Asfl = / \id\i\J DqDxD 2 
Jo 



(52) 



(53) 



Hence, for similar barriers we obtain the dependence 
J oc I? 3 / 2 rather than J oc D (as it would be the case 
for independent barriers). The latter dependence would 
follow from the calculation based on Zaitsev boundary 
conditions for the Eilenberger propagators. We observe, 
therefore, that quantum interference effects decrease the 
Josephson current in systems with three insulating bar- 
riers. This is essentially quantum effect which cannot 
be recovered from Zaitsev boundary conditions even in 
the multichannel limit. This effect has exactly the same 
origin as a quantum suppression of the average normal 
transmission (T) due to localization effects. Further lim- 
iting expressions for short junctions can be directly re- 
covered from eq. fl44|). 

We also note that the current-phase relation (|5^) devi- 
ates from a pure sinusoidal dependence even though all 
three transmissions are small -Do, 1,2 « 1 At T = the 
critical Josephson current is reached at \ — 1-7 which 
is slightly higher than 7r/2. Although this deviation is 
quantitatively not very significant, it is nevertheless im- 
portant as yet one more indication of quantum interfer- 
ence of electrons inside the junction. 

Finally, let us turn to the limit of long junctions 4^ » 
£o ■ We again restrict ourselves to the case of low trans- 
parent interfaces. At high temperatures T 3> Vp /2itd\^ 
from eq. (|l6|) , (|37|) we get 



J 



eTk 2 F A 2 sin x 

7T A 2 + 7T 2 T 2 



dmiD§D\Die 



(54) 



where d = d\ + d 2 and £(T) = vf/(2ttT). In this case 
the anomalous Green function strongly decays deep in the 
normal layer. Hence, interference effects are not impor- 
tant and the interfaces can be considered as independent 
from each other. In the opposite limit T — > (more pre- 
cisely T <C DvF/d), however, interference effects become 
important, and the current becomes proportional to D 5 / 2 
rather than to D 3 . Explicitly, at T — > we get 



J 



ek F VF sinx 
16n 2 y/did 2 Jo 



dnn 2 D 1 D 2 \/ D In D 



(55) 



This expression is valid with the logarithmic accuracy 
and no distinction between lnZ?o> InZ?! or ln_D 2 should 
be made. We see that, in contrast to short junctions, 
in the limit of thick normal layers interference effects in- 
crease the Josephson current as compared to the case of 



independent barriers. The result (pq), as well as one of 
eqs. ( f52| ) ( (53| ) cannot be obtained from the Eilenberger 
approach supplemented by Zaitsev boundary conditions. 



V. DISCUSSION AND CONCLUSIONS 

Let us summarize our key results and observations. 

In the present work we considered an interplay between 
the proximity effect and quantum interference of elec- 
trons in hybrid structures composed of normal metallic 
layers and superconductors. Quantum interference ef- 
fects occur between electrons scattered at different metal- 
lic interfaces or other potential barriers and can strongly 
influence the supercurrent across the system. 

The standard quasiclassical approach which describes 
scattering at interfaces by means of the nonlinear bound- 
ary conditions |0| for energy-integrated Eilenberger 
propagators - while very efficient in numerous other sit- 
uations - is in general not suitable for the problem in 
question. Because of this reason we made use of an al- 
ternative quasiclassical approach which allows to inves- 
tigate superconducting systems with more than one po- 
tential barrier and fully accounts for the interference ef- 
fects. Within this approach scattering at boundaries is 
described with the aid of linear boundary conditions for 
quasiclassical amplitude functions. Electron propagation 
between boundaries is described by linear quasiclassical 
equations. Our approach is technically not equivalent 
to one based on the Bogolyubov-de Gennes equations. 
In particular, our method allows to explicitly construct 
two-point Green functions of the system and bypass such 
intermediate steps as finding an exact energy spectrum of 
the system with subsequent summation over the energy 
eigenvalues inevitable within the Bogolyubov-De Gennes 
approach. On the other hand, if needed, the full infor- 
mation about the energy bound states can easily be re- 
covered within our technique by finding the poles of the 
Green functions in the Matsubara frequency plane. 

Within our method we evaluated the dc Josephson 
current in SNS junctions containing two and three in- 
sulating barriers with arbitrary transmissions, respec- 
tively SINI'S and SINI'NI" S junctions. For the sys- 
tem with two barriers and few conducting channels we 
found strong fluctuations of the Josephson critical cur- 
rent depending on the exact position of the resonant 
level inside the junction. For short junctions d -C £,oD 
at resonance the Josephson current does not depend on 
the barrier transmission D and is given by the stan- 
dard Kulik-OmePyanchuk formula |E| derived for ballis- 
tic weak links. In the limit of long SNS junctions d ^S> Co 
resonant effects may also lead to strong enhancement of 
the supercurrent, in this case at T — > and at resonance 
the Josephson current is proportional to D and not to 
D 2 as it would be in the absence of interference effects. 
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It is also interesting to observe that, while the above re- 
sults for few conducting channels cannot be obtained by 
means of the approach employing Zaitsev boundary con- 
ditions, in the many channel limit and for junctions with 
two barriers the latter approach does allow to recover cor- 
rect results. This is because the contributions sensitive 
to the scattering phase are effectively averaged out dur- 
ing summation over conducting channels or, which is the 
same, during averaging of the current over the directions 
of the Fermi velocity. 

Quantum interference effects turn out to be even more 
important in the proximity systems which contain three 
insulating barriers. In this case the quasiclassical ap- 
proach based on Zaitsev boundary conditions fails even 
in the limit of many conducting channels. In that limit 
the Josephson current is decreased for short junctions 
(J oc D z / 2 ) as compared to the case of independent 
barriers (J oc D). This effect is caused by destructive 
interference of electrons reflected from different barriers 
and indicates the tendency of the system towards local- 
ization. In contrast, for long SNS junctions with three 
barriers an interplay between quantum interference and 
proximity effect leads to enhancement of the Josephson 
current at T — > 0: We obtained the dependence J oc Z) 5 / 2 
instead of J oc D 3 for independent barriers. We also dis- 
cuss some further concrete results which turn out to be 
quite sensitive to the details of the model. 



Finally, we note that in a very recent publication ]2S| ] 
Ozana and Shelankov analyzed the applicability of the 
quasiclassical technique for the case of superconducting 
sandwiches with several insulating barriers. For such sys- 
tems they also arrived at the conclusion that in the many 
channel limit the standard quasiclassical scheme based on 
the Eilenberger equations and Zaitsev boundary condi- 
tions effectively breaks down in the presence of more than 
two reflecting interfaces. These authors argued that in 
such cases this scheme disregards certain classes of inter- 
fering quasiclassical paths. This conclusion p9| is similar 
to one reached in the present paper for SNS structures. 
We would like to point out, however, that from our point 
of view the failure of the above scheme is not so much 
due to the quasiclassical approximation and/or normal- 
ization conditions employed within the Eilenberger for- 
malism. The problem is rather in the boundary condi- 
tions jjJI which disregard interference effects which oc- 
cur in the structures with several interfaces/barriers with 
transmissions smaller than one. 

We would like to thank J.C. Cuevas, D.S. Golubev, 
A. A. Golubov, M. Eschrig, A. Shelankov, G. Schon and 
U. Ziilicke for discussions and useful remarks. The work 
is part of the CFN (Center for Functional Nanostruc- 
tures) which is supported by the DFG (German Science 
Foundation). We also acknowledge partial support of 
RFBR under Grant No. 00-02-16202. 



APPENDIX A: 



Let us consider an SINI'S system and assume that the normal metal layer is located at — d/2 < x < d/2. It is 
convenient to choose the coordinate x' within the normal layer, — d/2 < x' < d/2. Then a general solution of eq. ([!]) 
(decaying at x — > — oo) in the left superconductor reads 

( ?+ M ) = ( -iA- 1 ) e-'^m + ( J /2j ) e^e-**W). (Al) 
Here k = Cl n /v x . The solution in the right superconductor can found analogously. We get 

( ) = ( feA ) + ( -ie- V 1 ) (A2) 

The above solutions contain four unknown functions f(x'), g(x'), n(x') and r(x'). These functions should be found 
by matching (Al), (A2) with the solution of eq. ([j]) in the normal layer. The latter has the form 

( G F+t\t') ) = ( -£ e[lkx ~ { " Q M]lX ~ X ' ) + e-^^e^h(x') ( I ) + e «^e- ik **j(x') ( J \ (A3) 

+e" nX / v *e ik * x k(x') K J + e-^ x / v *e- tk * x l(x') ^ ^ 

and contain four additional unknown functions h(x'), j(x'), k(x') and l{x'). The boundary conditions at two NS inter- 
faces provide eight equations which allow to uniquely determine all the above functions and, hence, the supercurrent 
in SINI'S junctions. These equations are specified below. 
Consider the left boundary. Making use of eq. Jl^) we find 
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h(x')e 2 ™* 
k(x')e~%% 



0L\ 



,-nd/2 



while eq. 



yields 



j(x')e «S 



(A4) 



- Kd /' 2 f{x'). (A5) 



Similarly, applying the boundary conditions at the right interface one gets 



ie lx / 2 7 



e- Kd t 2 n{x') = a 2 



e 2 - k(x') 



■ft 



e 2 ?* j(>') 
e Z(a;') 



and 



_i e -<x/2 7 -i 



-«d/2 



-— e" ifex:r e e 2 "" + e 2 «™ /i(x') 



e 2 ™* k(x') 



It is easy to see that the free terms in eqs. (A4)-(A7) are 



zi(x') 



ik x x'—(a! n x' /v x ) 



Z2(X') 



__^ — ik x x' J t-{uj n x'/v x ) 



(A6) 



(A7) 



(A8) 



Eqs. (A4)-(A7) can be trivially resolved and we arrive at the solutions for the functions h(x') and j(x') (which are 
only needed in order to evaluate the current) with the structure 



h(x') = Uizi(x') + Uiz 2 (x), j(x') = Vizi(x') + V 2 Z2(x r ), 
where Ux,2 and V\ 2 do not depend on x' . 

APPENDIX B: 



(A9) 



Consider the Eilenberger quasiclassical propagator 
which has a 2 x 2 matrix structure in the Nambu space 



g(p, R, Lu r , 



g(p,R,u n ) if(p,R,u n ) 
-if+(p,R,ui n ) -g(p,R,w n ) 



(Bl) 



Here p is the electron momentum on the Fermi surface 
and R is its coordinate. This quasiclassical propagator 
obeys the normalization condition g 2 = 1 or, equiva- 
lently, g 2 + ff + = 1. In addition, anomalous (/, / + ) 
and normal (g) Green functions obey important symme- 
try relations 

/+* (p, R, lu u ) = f(-p, R, «„) = f( Pl R, -w B ), (B2) 
g*(p,R,uj n ) = g(-p,R,v n ) = -g{p,R,-u> n ). 

The Eilenberger equations 0] can be written in a concise 
matrix form as 



ivpWg + tog — gui = 0, 



(B3) 



u> = (iuj n + -v F A)a z - A + — (g) + —(a z ga 



impurities, respectively. The angular brackets denote av- 
eraging over the Fermi surface. The matrix A(_R) incor- 
porates the superconducting order parameter A(U), d z 
is the Pauli matrix 



A = 



A 

-A* 



1 

-1 



(B4) 



The current density j(R) is defined as follows 

j(R) = -2nieTN(0) ^(v F (p)g(p, R, «„)) (B5) 

m 

N(0) is the density of states at the Fermi energy per one 
spin direction. 

The Eilenberger equations (B3) should be supple- 
mented by Zaitsev boundary conditions at metallic in- 
terfaces. These conditions have the form p| 



9a+ = ga- 



ga, 



g a [(l - D(p))(g s+ + g s -) 2 + (g s+ - g s „) 

= Dip) [g s +g s - - g s -g s +] ■ 



(B6) 



(B7) 



Here g s , a ip, R, u n ) = \gip,R,u) n ) ± g(p r ,R,u n )]/2, by 
where A stands for the vector-potential; r and r s are the p r we denote the reflected momentum. The subscripts ± 
elastic scattering time on nonmagnetic and paramagnetic in eqs. (B6), (B7) stand for the expressions on the right 
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(left) side of the interface, respectively. Finally, D(p) 
is the transparency coefficient of the boundary for the 
electron at the Fermi surface with the given direction of 
momentum. Eq. ( p36| ) results in the current conservation 
at the boundary. 

Consider the Josephson current in a clean SINI'S 
structure in the absence of the magnetic field. We as- 
sume that both NS interfaces are specularly reflecting 
and are perpendicular to the cc-axis. In this case the 
quasiclassical propagator depends on p x , x and u> n . Mak- 
ing use of the symmetry relations ( p32| ) the functions 
9s,a(px,x,u> n ) = \g(p x ,x,u> n )±g(-p x ,x,u> n )]/2 and sim- 
ilarly defined functions f s , a , ft a can be parametrized as 
follows 

g s = f s = b s2 - ib s3 , /+ = b s2 + ib s3 , (B8) 

9a = ibal, fa = b a 3 + ib a 2, fa = ~ b a3 + «&a2, 

The parameters b s , b a are real and constitute two three- 



dimensional vectors b. 



(a) 



(b 



sl(al)) u s2(a2): &s3(a3))- 



Combining the Eilenberger equations for g(±p x , x, ui n ), 
one easily finds 



db s db a 
— = Mxb a , — 

ax ax 



-Mxb a . 



(B9) 



Eqs. ( |B9D should be considered only for positive p x > 
and oj n > 0. The three-dimensional vector M in ([B9| ) 
is real and has the following components Mi = 2uj/v x , 
M 2 = (A+A*)/v x and A/ 3 = i(A-A*)/v x . Note that by 
introducing a complex vector z = b s +ib a one can rewrite 
eqs. (B9) as dz/dx — —iM x z. From the latter equa- 
tion we conclude that z 2 should be equal to a constant. 
From the normalization condition one finds z 2 = 1, or 



b s = e M+ + C + exp(-|M|ar), b a = -b s x e M+ , 



(B14) 



where vector ejv/+ is given by eq.(B13) with the changed 
sign of X- Using these equations, from eq. (Bll) one can 
establish the relation between b a and b s at the normal 
side of the interface near the left superconductor 



tib s x e 



M- 



Similarly, for the right boundary we get 
b a = t 2 eM+ x b s . 



(B15) 



(B16) 



With the aid of these conditions one can easily find the 
Josephson current in SINI'S junctions. What remains 
is to solve the Eilenberger equations ( |B9| ) in the normal 
metal. In the absence of external fields and impurities 
we have 



C 

b„ = | L + coshi + L_ sinh.T 
M+ cosh x + M- sinh x 



D 

b a = I M + sinh x + M_ cosh x 
-L + sinh x — L- cosh x 



(B17) 



(B18) 



where x = 2to n x/v x and C, D, L±, M± are constants de- 
termined from the normalization and boundary condi- 
tions. Finally, making use of eq. (B5) we arrive at the 
result 



b 2 s = 1 + b 2 a , b s b a = 0. 



(BIO) 



The boundary conditions (pa), (B7) take the form 



b a - = b a+ = b a , (Bll) 

[(6,+ - 6 S -) 2 + (1 - D)(b s+ + 6 s _) 2 ]b a = 2Db s+ x b s 

Assuming the pairing potential to be constant in the 
superconductors one easily recover the solution of the 
Eilenberger equations. For the left superconductor x < 
—d/2 we obtain 



b s = bm- + C— exp(|JVf |x), b a = b s x e 



M- 



(B12) 



Here C_ is an arbitrary vector perpendicular to M and 
ejv/- is the unit vector in the direction of M 



eu- = 



Acos( X /2) 



vW+^n y Asin( X /2) 
Analogously, for x > d/2 we have 



(B13) 
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